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Abstract: We present supersymmetric scenarios with gravitino LSP and stau NLSP in 
the case of a non-standard model of cosmology with the addition of a dark component in 
the pre-BBN era. In the context of the standard model of cosmology, gravitino LSP has 
drawn quite some attention as it is a good candidate for dark matter. It is produced in 
scattering processes during reheating after inflation and from the decay of the stau. With a 
long lifetime, the stau decays during Big Bang Nucleosynthesis. It is strongly constrained 
by the abundance of light elements but can however address the known "BBN lithium 
problem". It requires fairly massive staus mf > 1 TeV and puts an upper bound on the 
reheating temperature Tr ~ 10 7 GeV which does not satisfy the requirements for thermal 
leptogenesis. For the non-standard cosmological scenario, the reheating temperature 
bound can be strongly relaxed Tr S> 10 9 GeV and the lithium-7 problem solved with a 
stau typical mass of m,f ~ 600 — 700 GeV and down to ~ 400 GeV with a very important 
dark component that could enable possible production and detection at the LHC. 
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1. Introduction 



From a cosmological perspective, super symmetry is a promising theory beyond the Stan- 
dard Model. Assuming the conservation of R-parity, the lightest supersymmetric particle 
(LSP), such as the neutralino or the gravitino, can be a candidate for dark matter. The 
gravitino acquires mass through the super-Higgs mechanism in broken local supersymme- 
try and which spans from few eV to TeV depending on the scenario. The gravitino LSP 
has been studied intensively @, |, § |, |, § 0, |, g 0, [TT], [12], [E| [U], [ll| [l§ [17|, 0, @ |0| 
as a candidate for dark matter. 

The current dark matter density has been measured with great precision by the WMAP 



mission [21], the five-year data giving at 3<r 



n DM h 2 = 0.1099 ±0.0124 



(1.1) 



where h is the reduced Hubble parameter. This result puts strong constraints on the 
production of gravitinos in the early Universe. Its relic density must lie below the upper 
limit to avoid an overclosure of the Universe. The lower bound on the dark matter is less 
stringent as one can have multi-component dark matter with gravitinos and other species. 

There are two main gravitino production processes. On the one hand, the thermal 
production (TP) controlled by the post-inflation reheating temperature produces graviti- 
nos in scattering processes during reheating ||, ||, ||, ||, |{| |7|, 11' H U> H HI 



32, 33, 18, 2C]. To avoid overclosure of the Universe, one obtains upper limits on the re- 
heating temperature Tr < 10 6-8 GeV. The only experimental constraint on lower bounds 
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for the reheating temperature comes from Big Bang Nucleosynthesis that can occur only 
if Tr > 1 MeV. But the derived bound is somewhat below the required temperature 
for thermal leptogenesis Tr > 0.4 — 2 x 10 10 GeV [34]. On the other hand, gravitinos 
are produced from the decay of supersymmetric particles. Since gravitinos only interact 
gravitationally with the other particles, these interactions are suppressed by the Planck 
mass and are therefore very small (in the limit of our considered gravitino masses, for very 
light gravitinos, the goldstino component increases the interactions). While the Universe 
cools down, SUSY particles would decay preferably to the next to lightest supersymmet- 
ric particle (NLSP) which would then decay to the gravitino. This process is called the 
non-thermal production (NTP) of gravitinos 0, [2|, [3|, EL |5j, §. The contribution of both 
NTP and TP processes contributing to the total gravitino relic density was performed in 
0, 12, 19, ^(J with the addition of Big Bang Nucleosynthesis (BBN) constraints to solve 



the lithium problems. In standard cosmology, the Big Bang Nucleosynthesis (BBN) is 
a very predictive model yielding the abundance of light elements produced in the early 
Universe. These predictions are confronted to observations and are in great compatibility 
for deuterium, in quite good agreement for helium-4 but present strong discrepancies for 
lithium-7 and lithium-6. Solutions to the lithium problems are investigated on an astro- 
physics perspective or on the particle physics side by assuming the decay of relic particles 
during BBN. Solutions for the lithium problems were studied in specific models such as 
CMSSM or GMSB models with stau or neutralino NLSP decaying to the gravitino LSP 
lH HI PI,® II- The NLSP is long-lived due to the Planck mass suppressed coupling 
to the gravitino and can decay during the BBN at times above one second. In the CMSSM 
framework with gravitino LSP and stau NLSP, it was shown that it is possible to solve 
both lithium problems and obtain the dark matter relic density from gravitino thermal 
and non-thermal production |3(| [19[ pCfl . The downsides are the necessity for quite small 
reheating temperature and a fairly heavy spectrum. 



In this study we consider CMSSM scenarios with a stau NLSP, a gravitino LSP pro- 
duced from TP and NTP and we assume a non-standard cosmological history. We introduce 
a modified expansion rate due to the energy density of a dark component in the pre-BBN 
era following the model presented by Arbey and Mahmoudi fl37| |38| , The consequence of 
such a modification of the Hubble parameter is a less efficient production of gravitino in 
the thermal process and an early freeze-out of the NLSP leading to higher abundance of 
the NLSP. The constraints on the reheating temperature are therefore less stringent and 
allow higher values compatible with thermal leptogenesis. The higher abundance of NLSP 
allows to solve the lithium problems with lighter masses for the NLSP than in the stan- 
dard scenario. Consequently these supersymmetric spectra could be accessible to collider 
experiments. 

This paper is organised as follows, in section [2] we describe the model with a pre-BBN 
expansion rate. In sections ^ and EJ we present the calculation of the gravitino relic density 
through thermal and non-thermal production followed in section || by the resolution of the 
lithium problems from the decay of long-lived stau. 
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2. Pre-BBN modified expansion 



Observations of the early Universe are mostly reliable for temperatures up to T ~ 1 MeV 
which correspond to the beginning of the Big Bang Nucleosynthesis. From that moment the 
Universe cooled down and its evolution is described by the standard model of cosmology, 
the hot Big Bang model. It is constrained by many observations such as the abundance 
of light elements, the cosmic microwave background radiation or the formation of large 
scale structures. The Universe must be dominated by radiation between BBN and the 
recombination followed by an era dominated by matter. On the contrary the evolution 
of the Universe in the pre-BBN era is unclear. There are no strong constraints on its 
composition. This open question can be related to the nature of dark energy, responsible 
today for the accelerated expansion of the Universe which is still unknown and could play 



an important role in the pre-BBN era. This is indeed the case in dark fluid models [ 39 J 



or quintessence models [40|. These dark energy components contribute to the Hubble 
parameter leading to faster or slower expansion rates (negative effective energies can be 
obtained in extra-dimension models |4l|l). 



To model the effects of a modified expansion in the pre-BBN era, we will use the 



parametrisation given by Arbey and Mahmoudi 37] adding a new dark density to the radi- 



ation density depending on temperature and characterising a fluid in adiabatic expansion 

/ rp \ n D 

p D (T) = PD (T BBN ) (2.1) 

V-'BBN/ 

where no is a constant describing the density behaviour and the reference temperature 
Tbbn = 10 MeV is chosen higher than the beginning of Big Bang Nucleosynthesis and 
neutrino freeze-out in order to assure a radiation dominated era for BBN. The equation of 
state \swd = Pd/pd where Pjj is the pressure of the fluid, tid and wd are related by no = 
3(wd + 1) (more details can be found in |4(|). The value rip = 3 (wd = 0) corresponds to a 
matter density, no = 4 (wd = 1/3) to a radiation density and njj = 6 (wd = 1) corresponds 
to quintessential kination (quintessence with kinetic energy dominating potential energy 
|4^| ). Higher values of njj can be reached in models with a decaying scalar field. 
As we require a radiation dominated era for BBN, we introduce the parameter 

Pd(Tbbn) , (r>r> , 

KD = ~Fr V < 1 t 2 - 2 ) 

/WUbbn) 

where /3 ra d is the radiation density 

p r ad(T) = <7efr(T)^r 4 (2.3) 

g e fi is the effective number of radiation degrees of freedom. Higher values of K£> increases the 
length of the dark component domination era. For kd — > 1, radiation and dark component 
are co-dominant at Tbbn leading to perturbations in the production of light elements. 

For the two parameters controlling this toy model of dark component, we will restrict 
our study to values < K£> 1 and 4 < n£> < 8 as in [37|. 
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The Friedmann equation reads 



H 2 = —(p B + PD ) (2.4) 

where the energy of the background ps = Prad + Pm was dominated by radiation in the 
past and is now dominated by matter. For our purpose ps — Prad- Furthermore once BBN 
starts, the dark component becomes negligible and will have little effect on the expansion 
rate from that time on. The addition of the dark component in Eq. ( |2.4[ ) will have important 
consequences for the calculation of the relic density of dark matter. 



3. Thermal production of gravitino 

During the reheating era after inflation, scattering processes producing one gravitino can 
become very efficient. Bolz et al. have considered a consistent thermal field theory 
approach for supersymmetric strong interactions. Pradler and Steffen have taken into 
account the full Standard Model gauge group SU(3) C x SU{2)l x U{1)y using the same 
procedure [31, [33 1. Another approach |3^] leads up to a factor 2 difference with the previous 



study. Although, this approach is more complete, we will follow the procedure used by 
Pradler and Steffen for implementation purposes. 

The Boltzmann equation describing the evolution of gravitino abundance reads 



dn 3/2 



3Hn 3/2 = C 3/2 (3.1) 



dt 

with 723/2 the gravitino number density and C3/2 the collision term describing the produc- 



tion and destruction of gravitinos. The Hubble parameter Eq. (^4) can be rewritten 



where M p \ = 2.4 x 10 18 GeV is the reduced Planck mass. Using the definition Y3/2 = n 3 / 2 / s 
where s is the entropy and assuming that previous to reheating all gravitino density is 
diluted by inflation, we have Y 3 / 2 (Tr) = with Tr the reheating temperature. Integrating 
the Boltzmann equation between Tr and TbbNi an d assuming that all thermal production 
of gravitino is negligible at BBN times leads to 

fTBBN C 3/2 (T) 

h ^Ttmfyf <3 ' 3) 

The entropy density is 

2vr 2 

s = -rrh eS T 3 . (3.4) 
45 

Thermal equilibrium is assumed for all supersymmetric particles during reheating therefore 
5eff = h cS = 915/4. 
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Pradler and Steffen [31[ have calculated the collision terms including scattering pro- 
cesses producing one gravitino in strong and electroweak interactions 



C3I2{ ) = 5 1 ^T 9 " t ) ■ ( ' 



The sum is over the gauge groups with the calculated coefficients c a = {11,27,72} and 
k a = {1.266,1.313,1.271}. The gaugino masses and the gauge couplings are respectively 
M a and g a . 

For kd = 0, we recover the standard expression for the Hubble parameter. One obtains 
the result obtained by Pradler and Steffen by solving Eq. (3.3) with Tbbn <C Tr 

£{»£■ 

a=l \ 67Jl 3/2, 



- £ ( 1 + S-) ^> {t) (ra) < 3 - 6 ' 



with y a = {0.65,1.6,4.3} x 10 -12 . One can note the dependence on the reheating tem- 
perature, the gravitino mass and on the gaugino masses. Higher temperature and gaugino 
masses or smaller gravitino masses lead to a larger production of gravitinos. 

For the case kd 7^ 0, we take into account the effects of the dark component on the 
Hubble parameter. We get 

- ( Ml\ 9 , k a f TR dT 



Y 3/2 (T BBN ) = 2.08 x 10- 23 £ 1 + c a gl In ^ / 

3m 3/2/ 9a JT 



9a JTBBN 



rp \«D-4 

1 + k d(t^ 



1/2 

(3.7) 



If no = 4, the integral is independent of T and the yield of gravitino is 



WTbbn) = — — Ips(Tbbn) (3.8) 

1 + K D 

The result is simply proportional to 1psC?bbn) and for «£>—)■ 1, we get Y" 3 / 2 (2~bbn) — >■ 
^ps(2bbn)/2, and for kd = we recover the result of the standard scenario. The contri- 
bution of the dark component to the Hubble parameter tends to weaken the efficiency of 
gravitino production. 

When n_o > 4, integration of Eq. fl3.7f ) can again be expressed with Eq. ([T(]) multiplied 
by an hypergeometric function 

Y 3/2 (T B bs) = IpsObbn) x 2 F X \l/N, 1/2; 1 + 1/N; -k d (Jr^j j (3-9) 

with N = no — 4. Since kd > 0, the hypergeometric function 2 F% takes values in the range 
]0, 1[ so the gravitino production is always smaller than the standard production. It can 
also be directly seen from Eq. (|3.7l). The relic density of thermally produced gravitino is 

njf 2 h 2 = 2.742 x 10 8 Y 3/2 (T BW ) (3.10) 
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T R (GeV) 

Figure 1: Thermal production for different gravitino masses and {kd^d} at mi/2 = 1 TeV and 
TO 3/2 = 1 GeV. The black line is the standard case corresponding to the maximal production of 
gravitino, the black dashed line has njj — 4. Then full lines have nr> — 5, short-long dashed line 
n d = 6, long dashed lines rip — 7 and short dashed lines no = 8, red, blue and green correspond 
respectively to k d = {10~ 8 , 1CT 5 , 1CT 2 } 

In all cases, we find that the gravitino production is proportional to the standard result 
calculated by Pradler and Steffen multiplied by a coefficient between and 1. The addition 
of the dark component will induce the reduction of the gravitino thermal production. The 
production is smaller as the dark component dominates, that is if one takes large values 
of kd or ri£). This can be easily understood: the Hubble parameter is much larger than in 
the standard scenario meaning that one exits the reheating period much faster leaving less 
time to produce gravitinos. 

The thermal production of gravitino Fig. [I] shows the gravitino relic density with re- 
spect to the reheating temperature at m^j = 1 TeV and a gravitino mass m 3 / 2 = 1 GeV. 
We recall from Eq. (|3.6| ) that the thermal production is more efficient for lighter gravitino 
masses. The grey area corresponds to the WMAP value of dark matter relic density and 
puts constraints on the reheating temperature. To avoid an overproduction of dark matter 
in the standard cosmological scenario (black continuous line), the reheating temperature 
must lie below Tr ~ 10 6 ~ 7 GeV. The constraint becomes more severe for lighter grav- 
itino masses but can go up to Tr ~ 10 9 GeV for a gravitino mass of m.3/2 — 100 GeV 
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and mi/2 = 1 TeV. Higher values of m\/2 also increase the production of gravitinos. 
These limits on the reheating temperature are quite low for leptogenesis requirements 
(Tr ~ 2 x 10 9 GeV). One must either turn to other mechanisms for baryogenesis (for 
instance ^3|) or limit the production of gravitino. This can be achieved in our model by 
inducing a modification of the expansion rate. Models with different values for kd and 
rt£> are presented Fig. [l]. As a general conclusion, they all produce less gravitinos than 
the standard scenario, as expected, and therefore weaken the constraint on the reheating 
temperature. The black dashed line corresponds to models with = 4 as in Eq. ( |3.8|) in 
the limit case of kd — > 1. The gravitino relic density is proportional to the Pradler-Steffen 
solution with a factor between 0.5 and 1, so the constraints on the reheating temperature 
are still close to the standard case. The gravitino production will be reduced with respect 
to the amount of dark component. For no = 5 and kd = 10 -5 the reheating temperature 
limit is ~ 10 9 GeV and for higher values of ud the thermal contribution to the relic density 
is negligible, yielding no limit on the reheating temperature. It is now possible to satisfy 
the lower bound on the reheating temperature deriving from thermal leptogenesis. 

Gravitino can also be produced from other mechanisms that must be added to the 
suppressed thermal production in order to obtain the dark matter relic density. 



4. Non-thermal production of gravitino 

This process comes from the decay of supersymmetric particles into the gravitino produced 
in the early Universe before the beginning of BBN. Due to the Planck mass suppression 
in the gravitino couplings, SUSY particles decay to the NLSP which then decays to the 
gravitino. One can assume a standard calculation of the abundance of the NLSP as if it 
was stable by solving the Boltzmann equation ||44|| . Then one can assume that each NLSP 
will produce one gravitino and the gravitino relic density is related to the abundance of 
NLSP scaled by the mass ratio 

NTP^2 = ^3/^^ (41) 
1 "1NLSP 

The unknown abundance of NLSP needs to be calculated. Its computation is based on the 
solution of the Boltzmann equation 



dn 
~dt 



3Hn - (a eS v) {n 2 - < q ) (4.2) 



where n is the number density of all supersymmetric particles, n eq is their equilibrium 
density, and (o" e fff ) is the thermal average of the annihilation rate of the supersymmetric 
particles to the Standard Model particles. The numerical computation has been automated 
under the standard scenario in programs such as MicrOMEGAs [45] or DarkSUSY [46] with 
absence of entropy production and decay of relic particles. 

The Boltzmann equation Eq. ( [4.2; ) shows that the modification of the Hubble parame- 
ter with the introduction of the dark component will have strong impact on the freeze-out 
of the NLSP and its abundance. This has been studied in the case of the neutralino LSP 



1 42]. Since the Hubble parameter is larger with the dark component, the freeze-out occurs 
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earlier leaving a larger relic density because less annihilations were performed. Similar 
conclusions are drawn in the stau NLSP case. 

We use the usual ratio x = ttinlsp /T of the NLSP mass over temperature and com- 
bining Eq. ( pO] ) and Eq. ( pL^ ) it yields 

^ = -# ffy2m 9 NLSP f 1 + -^fU (a eff ^) (F 2 - F 2 ) (4.3) 

with 

5 y 2 = ^(i+^^) (4.4) 

The standard result is obtained with p£> —> 0. 

The abundance of the NLSP is then obtained by integrating Eq. ( [4. 3D between x = 
and x = tonlsp/^o where To = 2.726 K is the temperature of the Universe today. The 
yield is related to the relic density as in Eq. ( |3.10| ) 

Q NLSP h* = 2.742 x 10 8 (yg^) Y(T ) (4.5) 



Combining Eq. (4.1) and Eq. fl4.5|) , the non-thermal part of the gravitino relic density 
is proportional to the gravitino mass and to Y(Tq). The latter is numerically computed 
using a modified version of MicrOMEGAs 2.2 @ with SuSpect 2.41 J|7|. 

In the context of the CMSSM, the NLSP is mainly the lighter stau or the lightest 
neutralino. In this paper we will only focus on the stau case. First we will recall different 
results concerning the relic density of stau in the standard scenario without dark compo- 
nent. CMSSM depends on a small number of parameters, namely m^j the gaugino mass at 
the GUT scale, mo the scalar mass at the GUT scale, Aq the trilinear coupling at the GUT 
scale that we take Aq = 0, tan j3 ratio of the Higgs vacuum expectation values and the sign 
of fj, the Higgs mass parameter. We scan the parameter space with G [100, 6000] GeV, 
rriQ G [10, 2000] GeV, A$ = 0, fj, > and tan /3 = 10 or 50. Fig. || presents masses and 
abundances of the stau (and neutralino) in this model assuming a standard cosmological 
history. The green area noted WMAP indicates a relic density compatible with Eq. (1.1). 
It is mostly relevant to the neutralino as it could be a candidate for dark matter if the 
gravitino is heavier. For the stau NLSP Fig. || emphasizes with respect to Eq. ( [4.1| ) that 
for lower values of mi/2, NTP alone cannot produce dark matter and one must consider 
TP contribution. We have also added exclusion contours from LEP coming mainly from 
the lower bound on the chargino mass and at tan f3 = 50 a region where renormalization 
group equations do not give correct solutions for the SUSY masses. We have also indicated 
limits on the NLSP mass. Roughly, at mjrt S> mo the stau is the NLSP. Its mass increases 
with the gaugino mass. The stau mass can be analytically approximated by [48| at weak 
scale 

m 2 f = ml + 0.15m 2 /2 - 0.23 cos(2/3)Mf (4.6) 
The relic density the stau would have if it had not decayed to the gravitino reads 



n,h 2 = (2.2 - 4.4) x 10- 1 {-p^r) (4-7) 
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CMSSM tan/9 = 10 CMSSM tan/3 = 50 




Figure 2: Contours of ^nlsp^ 2 in the mo — mi/2 plane of the CMSSM for tan/3 = 10 (left) and 
tan/? = 50 (right). In green the region where riNLSP^ 2 is compatible with WMAP measurements. 
In light pink the region where ^nlsp^ 2 < 0.091, in pink 0.128 < JInlsp^ 2 < 0.5, in red 0.5 < 
^nlsp/i 2 < 1 and in dark red f^NLSpft 2 > 1- The grey area noted "LEP" is an exclusion contour 
coming from lower bounds on SUSY masses from LEP and "tachyonic" signals that there are no 
RGE solutions. The dashed lines are contours of the NLSP masses. 



The coefficient in Eq. fl4.7[) depends on the importance of slepton coannihilation. In the 
region my 2 3> mo, it is clear that the stau relic density increases with m^ 2 independently 
of ttt-o. Dependence on mo becomes important in the stau-neutralino coannihilation region. 
For a given value of mj /2, the abundance of stau increases due to a smaller cross-section 
from coannihilation processes. At tan/3 = 50, the parameter region with stau NLSP is 
wider and the relic density is smaller due a larger left-right mixing and enhanced stau- 
Higgs coupling. 

Adding the dark component, we increase the value of the Hubble parameter with a 
positive effective energy density. Therefore we expect the freeze-out to occur earlier since 
the condition to freeze-out is approximately H >T with the reaction rate T oc n {av) . If 
the freeze-out occurs earlier, less annihilation occurs and the relic density must be higher. 

As in the standard scenario, Fig. [3| reproduces the relic density of the NLSP for tan j3 = 
10 (left) and tan/3 = 50 (right) and the dark component parameters are set to kb = 10 -5 
and rip = 8. 

Comparing with Fig. ||, one can immediately assess the increase in relic density for the 
NLSP. For example, a relic density of Q^h 2 ~ 0.1 was achieved for m 1( / 2 — 1500 GeV at 
tan/3 = 10 in the standard cosmological scenario and at m^ 2 — 900 — 1000 GeV in this 
modified cosmological scenario. The same observations can be made at higher values of 
tan/3. The direct consequence is that one produces the same density of NLSP for lighter 
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masses. Similar effect can be obtained for the neutralino NLSP. Constraints on dark matter 



calculated in Eq. (11) impose either lower values of m^j to compensate the increase in 
Clfh 2 or smaller gravitino masses. For a given mass, the increase in the stau abundance 
will have important consequences to solve the lithium problems as described in the next 
section. Fig. || represents the evolution of the stau relic density with respect to kq and 
no for rrii/2 = 1000 GeV, mo = 100 GeV and tan/3 = 10. The relic density increases 
very quickly when the density of the dark component becomes important with kd and no- 
Lines indicate the gravitino relic density compatible with WMAP from NTP only. The 
dark component contribution must be smaller for increasing gravitino mass. 

4.1 Total gravitino relic density 

The total gravitino relic density is the sum of thermal and non-thermal production 

n 3/2 h 2 = njf 2 h 2 + n^ 2 p h 2 (4.8) 

Both contributions are affected by the presence of a dark component acting on the expan- 
sion rate in the pre-BBN era. While increases Qj£ reduce. Fig. [| shows the limit 
on the reheating temperature in the standard scenario and in the modified expansion rate 
with kd = 10 -5 and no = 5. The areas correspond to £l 3 / 2 h 2 compatible with WMAP 
constraints. For small values of vn\i 2 i the relic density is dominated by the thermal pro- 
duction. Higher gravitino masses allow higher values of reheating temperature. In the 
standard scenario with 7713/2 = 100 GeV, one can observe the beginning of a change in 
the slope at mj/j ~ 5 TeV when NTP becomes dominant and alone reproduces all dark 
matter putting severe limits on the reheating temperature to ensure negligible thermal 
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Figure 4: Relic density of a stau NLSP in the kjj — no plane for tan/3 = 10, mxji = 1000 GeV 
and mo = 100 GeV corresponding to a stau mass mf ~ 370 GeV. Lines indicate gravitino relic 
density from NTP compatible with WMAP bounds on dark matter relic density for gravitino masses 
to 3/2 = 10, 100 GeV 

production. In the non-standard scenario, bounds on Tr are much higher because of the 
suppression illustrated in Fig. The suppression increases with Tr and as the limit on 
the reheating temperature are higher for low m\/2 in the standard case, the bounds in the 
non-standard model are much higher at low values of m\/2- Suppression is less efficient 
at low reheating temperature leading to a convergence of both scenarios at large values of 
mi/2- Bounds on Tr in the non-standard cosmologal picture are easily in agreement with 
thermal leptogenesis requirements. 

5. Lithium problems 

Big Bang Nucleosynthesis describes the production of light elements during the early Uni- 
verse. The success of the predictions of the standard model of BBN (SBBN) compared to 
the observations was a strong argument in favour of the hot Big Bang scenario as it implies 
that the Universe has expanded and cooled down from a very hot state and also suggests 
that only a small component of matter is baryonic, the rest being mostly dark matter. 
The SBBN model depends on a unique parameter, the baryon-to-photon ratio measured 



11 



CMSSM tan = 10 m =80 GeV 




Figure 5: Constraints on reheating temperature for f2 3 / 2 /i 2 with both TP and NTP contribution 
compatible with WMAP measurements, for the standard case and for a dark component with 
Kd = 10~ 5 and ud = 5 with gravitino masses m 3 / 2 = 1> 10, 100 GeV 



by WMAP mi 



m = 10 10 — = 6.225 ±0.170 (5.1) 



and therefore makes simple prediction on the abundance of light elements. Comparison 
of predictions and observations are nicely compatible for deuterium and helium-4, and 
somewhat inconclusive for helium-3 (due to important uncertainties in the post-BBN evo- 
lution of this element). However important discrepancies arise for lithium-7 and also for 
lithium-6. The lithium-7 abundance is inferred from absorption lines in the atmosphere of 
low-met allicity halo or globular cluster stars. A spectacular result namely the Spite plateau 
pQ[| put forward that, for low-metallicity stars, the Li abundance depended linearly with 
a small slope on metallicity and with very little scatter indicating that depletion would 
not have been very effective. One can interpret the extrapolation at zero-metallicity as the 
primordial abundance of Li. Many observations [51, 52 1 give an abundance in a range of 
'Li/H = (1.1 — 1.5) x 10~ 10 . Some uncertainties on the estimation of the atmospheric tem- 
perature of the stars could lead to higher estimations 7 Li/H = (2.19 ±0.28) x 10~ 10 [ 53 1 or 
even higher values 7 Li/H = (2.34±0.32) x 10" 10 @. The SBBN model prediction is a fac- 
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tor 3-4 higher than the observations. The latest estimations are 7 Li/H = 5.24^q' 67 x 10 10 
|55H and (5.14±0.50) x 1CP 10 f5"6fl . The measurement of lithium-6 has also lead to a possible 



6 Li problem. SBBN prediction are of order 6 Li/H ~ 10 15 — 10 14 [57] but substantial 
amount of 6 Li has been measured in stars 6 Li/'Li ~ 0.05 [ p2|j , most measurements corre- 
sponding to 6 Li/H = (3 - 5) x 10~ 12 |||. This lithium-6 issue is still controversial and 
more data is required to assure that there is a substantial production of lithium-6 during 
BBN. Here we will assume the possible existence of large amount of lithium-6 conflicting 
with its standard abundance which should be very low because the only 6 Li production 
reaction in SBBN D(ce,7) 6 Li is a quadrupole transition. Most 6 Li is produced by cosmic 
ray nucleosynthesis from spallation or cosmic ray fusion, but the abundance of 6 Li in low- 
metallicity stars is too high by two orders of magnitude to be explained by standard cosmic 
ray nucleosynthesis. The lithium discrepancies could find their origin in numerous sources 
(for a review see (5^]). First for lithium-7, the nuclear network used for BBN calculations 
could have some uncertainties thus overestimating 7 Li and underestimating 6 Li abundances 
|6C. 31 1 . It seems very unlikely that the uncertainties could bridge completely the discrep- 
ancies. A second hypothesis comes from possible systematic errors in the determination of 
7 Li and 6 Li abundances. Another possibility is that the measured abundance is not primor- 
dial but has been reduced by stellar depletion, atmospheric 7 Li could be destroyed while 
transported to inner parts of the stars by nuclear burning. However both the presence of 
the Spite plateau with little scatter and the presence of large amount of 6 Li which is more 
fragile than 7 Li and should be even more depleted argue against a simple depletion mech- 
anism. Some models with turbulent mixing |62|| can account for 7 Li. Combining atomic 
diffusion and turbulent mixing [63] can deplete efficiently 7 Li bringing SBBN prediction 
and observations in satisfying agreement but is somewhat artificial and ad hoc. 

A solution for both lithium problems can be found by modifying SBBN. Relic parti- 



cles could decay during BBN producing Standard Model particles p4L 65, 35, p§|. These 



particles would induce non-thermal reactions that would not disturb much the abundance 
of the other elements but could produce 6 Li through reactions such as 3 H(a, n) 6 Li or 
3 He(a,p) 6 Li with 3 H and 3 He produced from 4 He spallation. Non-thermal reactions, so- 
lutions for the Li problem were also investigated. Early attempts of reduction of Li 
by photo-disintegration || were incompatible with lower limit on D/H or upper bound on 
3 He/D [66, 67], but injection of neutrons during BBN reduces the 7 Li abundance [35]. This 
scenario has been studied intensively as it was a possible solution to the lithium problems 

m i, i, i, 0, uj, n n, uj, in p. 



In |L9|, Fig. 2 and Fig. 3 give model independent requirements for the NLSP abun- 
dance, hadronic branching ratio and lifetime to solve lithium problems. For the case of 
a stau NLSP with a typical hadronic branching ratio I?had ~ 10 -3 , the abundance must 
be of order Qfh 2 ~ 0.5 corresponding to a stau masses of 0(1 TeV). To get lifetimes of 
r ~ 3000 sec, the gravitino mass must be of order O(100 GeV). These rough estimates 
were confirmed in precise calculations in CMSSM or GMSB models UJI but this requires a 
quite heavy SUSY spectrum rendering almost impossible the production of SUSY particles 
at the LHC. 

By using the non-standard cosmological scenario with the addition of the dark com- 
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ponent, it is possible to satisfy the previous requirements to solve the lithium problems 
with lighter stau. From the construction of the dark component model, we impose that 
radiation dominates during the BBN era. Therefore we do not expect strong variations 
of the abundance of light elements produced during BBN due to a modified expansion 
parameter. It is confirmed by Arbey & Mahmoudi [38| in Fig.l in which the production of 



4 He and D are close to SBBN values in the correct parameter range for up and ri£). The 
authors observe strong productions of these elements for small values of no and high values 
of kd meaning that the dark component dominates during BBN which conflicts with the 
requirement of radiation domination during this period. 

To solve the lithium problems, we combine the known effects of a long-lived relic parti- 
cle decaying during Big Bang Nucleosynthesis and the non-standard cosmological scenario 
with a modified expansion rate. The larger Hubble parameter increases the abundance of 
relic particles, in our case, the stau NLSP decaying to the gravitino. And because of the 
larger abundance of stau, it is possible to solve the lithium problems with a lighter stau 
NLSP and a lighter gravitino. 

5.1 Calculation procedure 

The computation of the production of light elements during BBN must include decaying 
relic particles and bound state formation between the charged stau and the light elements. 



These bound states have catalytic effects on BBN and were first stressed by Pospelov [S8| 
focusing on lithium-6 production becoming very efficient through the reaction ( 4 He — f) + 
D — > 6 Li + r. Many other processes were studied 6£, 15, |7(| [72| and full cross-sections 
were calculated by Kamimura et al. [73]. The BBN computation is based on the SBBN 
Kawano code and modified by K. Jedamzik p0| , 7C] to include all the latest improvements 
on decaying relic particles and bound states. The Jedamzik code requires a number of 
parameters: the hadronic and electromagnetic branching ratios for the decay of the NLSP 
and the energies in the cascades. The NLSP decays to the gravitino and various Standard 
Model particles such as photons, electrons and positrons (electromagnetic cascades) or 
neutron and protons (hadronic cascades). These particles interact with the light elements 
present in the plasma inducing reactions like photo-dissociation, spallation. . . 

The stau decay is mostly dominated by the two-body decay f — > tG therefore its 
lifetime is given by the decay width of this process r = h/T(f — > tG) 

,2 \ 4 



5 / m 



r(f -»tG) = t~ "" T 2 1 1- 1 (5-2) 

Electromagnetic cascades dominate the whole decay process induced mostly by the 
unstable lepton tau produced in two-body decay. Below r ~ 100 s, the mesons produced 
in the decay of the lepton tau induce charge-exchange reactions increasing the neutron-to- 



proton ratio and the helium-4 abundance ([28] and references therein). Such effects are 



not taken into account and therefore, only lifetime above r ~ 100 s are considered. Above 
a NLSP lifetime of 100 s, mesons mostly decay electromagnetically. 

Hadronic cascades require the production of a quark-antiquark pair that will then 
hadronize and produce protons and neutrons (for which the energy spectrum is computed 
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using PYTHIA [|74]]). Therefore 4-body diagrams are required for stau hadronic decay width 
calculation. The hadronic decay is therefore sub-leading with respect to the electromagnetic 
decay. 

The electromagnetic branching ratio and the energy (simply given from kinematics) 
are 2 

B em = 1 - B had ~ 1 and E em = a N ^ SP — (5.3) 

The lepton tau produced in the decay of its superpartner is unstable and decays producing 
lighter leptons and neutrinos. The part of energy taken by the neutrinos has no or little 
effect on the plasma, the coefficient a signs this loss of energy and following |J we have 
1/3 < a < 1. We take a = 1/2. 

The hadronic branching ratio is calculated following JHj 

5had(f -> rGqq; m°f) = F(r ^ (5 . 4) 

J- tot 

with the 4-body decay width 



m ,-m 3/2 - m dr(f -» rGqq) 
Y[t -> TGqq;m q q) = / (5.5) 

where m 9 q = 2 GeV is a cut on the invariant mass of the pair quark-antiquark corresponding 
to twice the mass of a nucleon below which hadronization does not produce any nucleon. 



The decay width is calculated using CalcHEP [75] as in [19]. Diagrams taken into account 
are given in Fig. ||. The Z boson width is taken as a Breit-Wigner function. Diagrams 
with T2, neutralinos x^ X31 X4 an d charginos xf are neglected assuming these particles 
are very massive. Diagrams with the exchange of a virtual neutralino or chargino are 
negligible for m x > l.lmf but are always included in the calculation. Contributions to 
the hadronic branching ratio from diagrams with a virtual W are one to two orders of 
magnitude smaller than the contribution of processes with a virtual Z because the stau is 
dominated by its right component. Processes with Higgses are neglected because of the 
masses in the propagator and the small couplings. In our calculations, we fully take into 
account the left-right mixing of f\ and gaugino mixing of x?- 

The hadronic energy calculation is done following again [14] using a mean value based 



on a convolution of the branching ratio with the quark-antiquark invariant mass 

^had - — -~— / am q q m qq (5.b) 

r(r rGqq) J mqq dm qq 



It was mentioned in IS] that this estimation of the hadronic energy had some problems. 
First the energy is estimated in the rest frame of the quark-antiquark pair and not in the 
rest frame of the stau. This calculation does not allow to estimate the boost needed to 
change from one frame to another. Also the fact that this is a mean energy disregards 
completely the fact that one nucleon does not have the same impact on the plasma as 
ten nucleons with a total energy equal to the one nucleon. A more realistic generation 
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Figure 6: Diagrams calculated in the 4-body decay of stau to gravitino 



of nucleons and their energy in the rest-frame of the stau was used with CalcHEP and 
PYTHIA yielding an estimation of the uncertainties on the abundance of light elements. 
This study resulted in 10-20% uncertainty for deuterium, 30-40% for lithium-7 and 20% 
for lithium-6. In this paper we will use Eq. (5.6) keeping these uncertainties in mind. 



5.2 BBN results 

BBN calculation can now be performed taking into account a great number of effects: 
SBBN reactions, non-thermal reactions from the stau decay, bound state formation and 
the modified expansion rate. Abundance of deuterium, helium-3, helium-4, lithium-6 and 
lithium-7 are estimated and confronted to observational constraints. These constraints are 
taken on conservative basis to include observational uncertainties, post-BBN evolution and 
theoretical uncertainties in some reaction rates and hadronic energy calculation. We take 
the following BBN constraints : 

1.2 x 10~ 5 < D/H < 5.3 x 10~ 5 (5.7) 
3 He/D < 1.72 (5.8) 
Y p < 0.258 (5.9) 

8.5 x 10~ n < 7 Li/H < 2.5 x 10~ 10 (5.10) 

0.015 < 6 Li/ 7 Li < 0.66 (5.11) 

The last two constraints are the abundances satisfying lithium-7 and lithium-6 abundances 
inferred by observations. In the following figures, SBBN values for lithium isotopes will be 
taken as 

2.5 x 10~ 10 < 7 Li/H (5.12) 
6 Li/ 7 Li < 0.015 (5.13) 

Fig. |7| presents the light element abundances in the standard cosmological scenario 
with the CMSSM model and mo = 80 GeV and tan /3 = 10. Gaugino mass is chosen 
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m =80GeV tan|3=10 




Figure 7: Big Bang Nucleosynthesis results in CMSSM with mo = 80 GeV and tan/3 = 10. 
The light blue region produces light elements abundances compatible with SBBN, masses solving 
the lithium-7 discrepancy (yellow), lithium-6 (green) and both problems (red). The hatched area 
labelled NTP indicates a gravitino relic density from non-thermal production alone compatible with 
measurements from WMAP and the other hatched correspond to the sum of TP+NTP for given 
reheating temperature. A dashed line also indicates the 100 s NLSP lifetime limit. 

such that the stau is the NLSP so mi/2 = [400 — 6000] GeV and the gravitino mass 
m 3/2 = [10 _1 — 700] GeV. Strictly speaking, in the CMSSM the gravitino mass is of 
the order of the electroweak scale and lighter masses are more typical of gauge mediated 
supersymmetry breaking scenarios. It was shown in [|TJ]] that the choice of mediation 
mechanism did not change much the phenomenology of BBN therefore we use the CMSSM 
model and assume that the gravitino mass is a free parameter. The limits from deuterium, 
6 Li and 7 Li Fig. [7] delineate different regions. Light blue corresponds to values of mu2 and 
771,3/2 satisfying the SBBN values. Yellow indicates masses solving the 7 Li problem, green 
region satisfies the 6 Li corresponding to stau lifetime of 1600 — 6300 sec for an upper limit 
of 6 Li/ 7 Li < 0.66. For a more stringent limit of 6 Li/ 7 Li < 0.15, the stau lifetime does not 
exceed 4 x 10 3 sec. Above the 6 Li/ 7 Li > 0.66 limit, lithium-6 is overly produced because of 
bound state catalytic effects mainly the ( 4 He — r) + D — > 6 Li + f reaction. The red region 
solves both problems simultaneously for m 1 / 2 = [3 — 5] TeV corresponding to stau masses 
rrif = [1 — 1.8] TeV and gravitino masses = 60 — 120 GeV. Note that for such heavy 
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m„=80 GeV k d =1 s n D =5 m =80 GeV k d =1 s n D =8 




Figure 8: Results for Big Bang Nucleosynthesis in the modified expansion rate scenario with 
Kd = 10~ 5 and nu = 5 (left) or no — 8 (right). Colours and hatched areas defined as in Fig. |7[ 

spectrum, the gluino has a mass of almost 10 TeV which makes impossible the production 
of such particle in the LHC. For higher values of tan (3 would reduce the relic density of 
stau and therefore would require much more massive sleptons to produce the same impact 
of BBN. 

The hatched zone labelled NTP corresponds to the gravitino relic density coming only 
from non-thermal production compatible with WMAP constraints for dark matter. If one 
adds the contribution from thermal production, different regions satisfy the WMAP values 
depending on the reheating temperature. For Tr < 10 GeV, it is possible to satisfy to 
SBBN and dark matter constraints. In regions solving lithium problems, all dark matter 
is reproduced by gravitino relic density for Tr ~ 10 7 GeV. For smaller values of the 
reheating temperature, gravitino relic density lies below dark matter relic density and one 
must assume a mixed dark matter scenario with the gravitino and another candidate. At 
reheating temperatures Tr > 10 s GeV, constraints on gaugino mass and gravitino mass 
which must be heavier to limit the effect of thermal production do not permit to solve BBN 
constraints and are therefore excluded. 

Similar results are presented Fig. | in a non-standard scenarios with the dark com- 
ponent. We use the same CMSSM parameters and the dark component are taken to be 
kd = 10~ 5 and tid = 5 (left) and n D = 8 (right). For n D = 5, the effect of the dark 
component is strong enough to suppress the thermal production and relax the reheating 
temperature upper bound as mentioned before. It is now possible to solve both the lithium 
and dark matter problems with a reheating temperature up to Tr ~ 10 11 GeV, 4 orders of 
magnitude above the limit in the standard indicated Fig. |5| (for a gravitino mass 

m 3/2 = 1 GeV). For this scenario, effects on BBN is similar to what was found in the stan- 
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dard scenario. The dark component induces an early freeze-out of the stau and therefore 
increases its abundance but for these parameters, kd = 1CP 5 and no = 5, the variation on 
the stau abundance differs of less than a percent and the BBN results are unchanged. 

The figure on the right has an index more important no = 8 leading to an energy den- 
sity of dark component much more important. The consequence on the thermal production 
is a very strong suppression illustrated Fig. || (5 orders of magnitude at Tr = 10 5 GeV and 
the suppression increases with Tr). Only very little abundance of gravitino is formed and 
can be neglected with respect to the non-thermal production. There are no constraints on 
the reheating temperature any more. And one can choose any reheating temperature to 
be compatible with for instance thermal leptogenesis. 

On the other hand stau production is very efficient, enhanced by two orders of mag- 
nitude (see Fig. Q), because of the early freeze-out. The gravitino relic density from the 
non-thermal production is increased by the same factor. For BBN, the production of light 
elements changes strongly because much more unstable relic particles are present in the 
plasma. Upper bound on deuterium and lower bound on lithium-7 are reached for smaller 
values of m\j2- Bound state effects are also much more important since the abundance of 
stau is higher and limits on lifetime are much more severe: r < 2500 sec for 6 Li/ 7 Li < 0.66 
or r < 2000 sec for 6 Li/ 7 Li < 0.15. Abundances compatible with observations and the 
resolution of the lithium problems can be obtained with lighter stau than in the standard 
scenario rrif ~ 600 — 700 GeV and gravitino masses m 3 / 2 ~ 20 GeV. In this specific case of 
parameter choice for the dark component, it is not possible to solve both lithium problems 
and reproduce dark matter relic from gravitino alone. It is however possible to solve the 
lithium-7 problem and satisfy dark matter. In the region solving both lithium problems, 
the thermal production is negligible and the non-thermal production contributes to 10-20% 
of the dark matter relic density. One must keep in mind that dark matter could be mixed 
with another component. 

6. Conclusion 

Super symmetric models with gravitino LSP and stau NLSP have interesting applications 
for two cosmological issues, namely dark matter and Big Bang Nucleosynthesis. The grav- 
itino production from scattering processes during reheating and the decay of the NLSP 
can account for all dark matter in the Universe. And moreover, the long-lived stau decay 
brings solutions to both lithium problems. This nice set-up requires a very heavy spectrum 
and a quite low reheating temperature. 

In this paper, we have supposed that the expansion rate of the Universe previous to 
BBN differs from the standard picture because of the presence of a dark component. The 
consequences are various. On the one hand the thermal production of gravitino can be 
strongly suppressed relaxing the constraints on the upper bound for the reheating temper- 
ature. On the other hand stau freeze-out occurs earlier leading to a higher abundance for 
the slepton and a more important gravitino non-thermal production. And because of the 
increase in the stau abundance, equivalent impact on BBN as in the standard case can be 
reached with a lighter stau. 
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The very heavy spectra of the standard scenario is completely out of reach of the LHC, 
but the non-standard cosmology yield a lighter spectrum that may be investigated at the 
LHC. We have illustrated the case = 10~ 5 and rip = 8 where mf ~ 600 — 700 GeV 
which could be still quite hard to detect since squark masses are above 2 TeV. If kd = 
10~ 8 , mf ~ 400 GeV, supersymmetric events could be produced at the LHC running 
at 14 TeV. Detection was studied for instance in [76, 77]. Every supersymmetric event 



would induce a cascade ending with the lightest stau which will be stable on the scale 
of the detector (lifetime r ~ 10 3 sec) and would be identified as a heavy muon with low 
velocities and distinctive time-of-flight and energy-loss. The signature is quite clean to 
be distinguished from background and detection could give a measurement of the stau 
mass with its consequences for our dark matter and BBN scenario, illustrating a nice 
complementarity between cosmology and collider physics. 
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